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It has been suggested in Hoyng (2009) that dynamo action can be analysed by expansion of the magnetic field into dynamo modes 
and statistical evaluation of the mode coefficients. We here validate this method by analysing a numerical geodynamo model 
and comparing the numerically derived mean mode coefficients with the theoretical predictions. The model belongs to the class 
of kinematically stable dynamos with a dominating axisymmetric, dipolar and non-periodic fundamental dynamo mode. The 
analysis requires a number of steps: the computation of the so-called dynamo coefficients, the derivation of the temporally and 
azimuthally averaged dynamo eigenmodes and the decomposition of the magnetic field of the numerical geodynamo model into 
the eigenmodes. For the determination of the theoretical mode excitation levels the turbulent velocity field needs to be projected 
on the dynamo eigenmodes. We compare the theoretically and numerically derived mean mode coefficients and find reasonably 
good agreement for most of the modes. Some deviation might be attributable to the approximation involved in the theory. Since the 
dynamo eigenmodes are not self-adjoint, a spectral interpretation of the eigenmodes is not possible. 
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1. Introduction 

The origin of the magnetic field of the Earth is gener- 
ally understood in terms of dynamo action in the liquid 
outer core. Helical convection and shear flows generate 
and maintain the magnetic field against resistive decay. 
Hydromagnetic simulations confirmed this picture and 
many groups have published num erical geodynamo mod 
els in the las t 15 years (se e, e.g., iGlatzmaier and Robert: 



1995 



1997 
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Kuang and Bloxham, 1997; Kagevama and Sat- 
Dormv et al., 1998; Roberts and Glatzmaier, 2000 
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Christensen and Wichtl [2007). Even though the computational 
resources d o not permit the paramet e rs of the simulations to be 
Earth -like (IChristensen and Aubern. l2006t IChristensen et all 
2010l) . the magnetic field of the simulations reproduces many 
features of the observed geomagnetic field. They have also 
confirmed the idea of columnar convection, i.e. that the flow in 
the core is oft en organised in vortices spiralling parallel to the 



rotation axis (iBusse , 



1975). 



The simulations opened up the possibility for detailed di- 
agnostics of dynamo action. Parameter studies have been 
undertaken to classify the properties of the simulations and 
to find which parameter combi nations produce plane t ary or 
earth-like dynamos (see, e.g.. iKutzner and Christense nl 12001 
Christensen and AubertL l2006t IChristensen et all I2010I) . The 



mechanisms that generate a nd sustain the magne t ic fiel d have 
also been under scrutiny. Kag evama and Satol ( 1997 ) iden- 



tified elements of the aft mechanism of mean field theory, 
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while Olson et al. (1999) concluded that in their simulations 
the generation of poloidal and toroidal field resembles the a 2 
dynamo scenario. Another direction of research has been to 
try and id entify the mechan i sms that induce polarity reversals 
(see, e.g., Sarson and JonesUl999l: IGlatzmaier and Coe , 2007 : 
Dris coll and QlsonL 120091 : ISchrinner et all l2010al) . Reversals 
are sometimes initiated by a patch of reversed flux appearing 
deep in the core, which subsequently grows and spreads by 
advection over a larger region and eventually dominates over 
the old polarity fie ld. But no clear picture has emerged yet. 
Aubert et al.l y008) have developed a new visualisation tech- 
nique (dynamical magnetic field line imaging) and they find 
that reversals and excursions of the dipole axis are associated 
with what they refer to as magnetic upwellings. 

Another benefit of numerical dynamo models is that they per- 
mit to measure from the simulati ons the tensors an and fajk 



utilised in mean field theory, and Schrinner et al. (2007) have 
computed these m ean-field tenso rs with the help of the test field 
method. Recently. iHovng (2009) proposed a decomposition of 
the magnetic field in dynamo eigenmodes and computed the 
statistical properties of the multipole coefficients. The purpose 
of the present paper is to validate this expansion technique and 
to illustrate some of the capabilities and problems. 

To this end we decom pose the magnetic fi eld of the numeri- 
cal geodynamo model of lOlson et all d 19991) in dynamo modes 
and compare the statistical properties of the numerically derived 
mode coefficients with those derived by the theory of |Hovng 
(2009). The various steps involved are surveyed in Fig.Q] The 
numerical dynamo model is described in Section [2] The mode 
decomposition involves various steps presented in Section[3] 
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Figure 1: Survey of the computations reported in this paper. Starting point is the numerical dynamo model described in Sectionff] The top line illustrates how the 
field B is decomposed in dynamo modes. First we obtain the dynamo coefficients, then the dynamo eigenmodes and their adjoints. End product are the correlation 
coefficients {a'a k *). This computational train is treated in SectionfJ] The correlation coefficients may also be computed theoretically, which requires an additional 
computation outlined in the bottom line and described in Section|4] Measured and theoretical correlation coefficients are compared in Section lT2l 



In Section [4] we recapitulate the theoretical mode excitation 
levels and in Section [5] we compare the numerical results with 
the theoretical predictions. Finally, in Section|6l we summarise 
and discuss our results. 



grow exponentially if it is kinematically advanced with the 
time-dependent, saturated, velocity field taken from the self- 
consistent calculation. 



2. Numerical model 

We consider a Boussinesq fluid with electrical conductiv- 
ity a in a rotating spherical shell V and solve the momen- 
tum equation, the induction equation and a temperature equa- 



mperature e 

tion as described in more detail by Olso n et al 1 Jl999h . No- 
slip mechanical boundary conditions are applied and the mag- 
netic field continues outside the fluid shell as a potential field. 
Convection is driven by a temperature difference AT between 
the inner and outer spherical boundaries. The equations are 
governed by four dimensionless parameters, these are the Ek- 
man number E = v/QL 2 , the (modified) Rayleigh number 
Ra = aTgoATL/vSl, the Prandtl number Pr = v/k and 
the magnetic Prandtl number Pm = vjr\. Here, v denotes 
the kinematic viscosity, 57 the angular rotation rate, L the shell 
thickness, olt the thermal expansion coefficient, go the gravita- 
tional acceleration at the outer boundary, and n is the thermal 
and r) = 1 //z<7 the magnetic diffusivity with the magnetic per- 
meability fj, and conductivity a. 



The model under consideration is taken from iOlson et al 
1999b with E = 1(T 4 , Ra = 334, Pm = 2 and Pr = 1. The 



resulting convection pattern is strongly columnar and symmet- 
ric with respect to the equatorial plane, see Fig. [2] The magnetic 
field is largely dipolar, strictly antisymmetric with respect to the 
equatorial plane and approximately axisymmetric. The energy 
in the axisymmetric component of the field is about 40% of the 
total magnetic energy. 

The magnetic energy density exceeds the kinetic energy den- 
sity, as illustrated in Fig. [2] The magnetic Reynolds number 
Rm = u x , m . s ,L/r) is approximately 88. Here U T . m . s , denotes 
the r.m.s. magnitude of the velocity. Although kinetic and 
magnetic energy densities vary chaotically in time, the dipole 
axis of the magnetic field is stable and reversals do not oc- 
cur. This example belon gs to the class of kine matically sta- 
ble dynamos identified bv lSchrinner et al. (2010a), their model 
2, which means that there exists no magnetic field that will 
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Figure 2: Top: The radial velocity V r in the corotating frame of the dynamo 
model under consideration. The snapshot has been taken at r = 0.62 ro, where 
r u is the outer shell radius. The velocity component is normalised to its abso- 
lute maximum. Thus, the greyscale coding runs from —1, white, to +1, black, 
and the contour lines correspond to ±0.1, ±0.3, ±0.5, ±0.7, ±0.9. Bottom: 
Magnetic (upper curve) and kinetic energy densities of the model versus mag- 
netic diffusion time L 2 /r/. 



3. Mode decomposition 

In principle the magnetic field of a dynamo may be decom- 
posed in any suitable set of base functions. We could take, for 
example, the decay modes of a homogeneous sphere or spheri- 
cal shell. However, there are reasons to prefer another set: the 
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eigenfunctions of the dynamo operato r averaged ove r time and 



see 



Hoyng (2009) and Sec- 



azimuth of the numerical model, 
tion [5] The dynamo operator describes the average dynamo 
action of a given flow. For statistically steady and kinemati- 
cally stable dynamos, to which we restrict ourselves here, the 
actual magnetic field of the numerical model is expected to be 
well represented by the fundamental eigenmode and a few over- 
tones. Since the magnetic field of the numerical dynamo model 
is largely axisymmetric, we restrict the decomposition to the 
axisymmetric dynamo modes. Likewise we ignore symmetric 
modes since the field of the model is strictly antisymmetric with 
respect to the equatorial plane. 

A mode decomposition then requires three steps. The first 
is that we compute the dynamo coefficients ajj and ftijk of the 
model. Then we determine the eigenfunctions and eigenvalues 
of the dynamo equation, and finally we decompose the field B 
of the numerical model in these eigenfunctions. 

3.1. Properties and notations 

Before we engage in details we summarise some of the nota- 
tions and properties that we employ: 

- The magnetic field £>, vector potential A, current J and the 
velocity V in the dynamo region V are written in calligraphic 
font; B = V x A, etc. Furthermore, we make the usual split 
in mean flow v and turbulent flow u: V = v + u. 

- B is the mean field (B), and B l are the eigenfunctions of the 
dynamo equation of the numerical model, briefly referred to 
as the dynamo modes. We use upper indices to enumerate 
modes and lower indices for vector components. 

- The mode decomposition makes use of the adjoint eigenfunc- 
tions, indicated by a hat e.g. B h is the adjoint of B k . 

- The adjoint operation and V commute. The basic entities are 
the field and its adjoint. The adjoint current J k is defined as 

V x B k , likewise if B = V X (• • •), then (•••) = A. 

- We absorb the factor 4ir/c in the definition of all currents. So 
J = V x B, J k = V x B k , etc. 

- b l are the decay modes of the spherical shell V; the decay 
modes are self-adjoint, so we may adopt b = b % *\f — j" — 

V x 6", and we shall assume that they are orthonormal: 



b k * ■ b l d 3 r = S kt 



(1) 



V+E 



Here E is the vacuum exterior to V (not to be confused with 
the Ekman number). 

3.2. Dynamo coefficients 

The action of the small-scale turbulent flow u on the mean 
magnetic field is parametrised by the so-called dynamo coeffi- 
cients, a. and (3. They are, in general, tensors of second and 
third rank, respectively. The various coefficients describe, e.g., 
the anisotropic dynamo effects of helical flows (a-effect), the 
transport of the mean magnetic field and the anisotropic turbu- 
lent magnetic diffusivity. 

The dynamo coefficients are computed with the help of th e 




Figure 3: Two components of the dynamo coefficients. Left: a. vv in units of 
ri/L. On the right: D vv = rj + \fi^ r e — ^(3 v o r m units of t). For each 
plot the grey scale is separately adjusted to its maximum modulus with white 
as negative and black as positive. The contour lines correspond to ±0.1, ±0.3, 
±0.5, ±0.7, and ±0.9 of the maximum modulus. 



Table 1: Eigenvalues of the first 13 antisymmetric and axisymmetric dynamo 
mode s in units of r//L 2 , num bered according to decreasing 5ft A, and computed 
as in Schri nner et al] |2010b), with n max = 16, /max — 32. For a complex 
pair of modes, the larger mode number is taken to correspond to < 0. 



mode number i 


eigenvalue Ai 





(-3.875, 0.00) 


1,2 


(-34.83, ±10.31) 


3 


(-40.63, 0.00) 


4 


(-42.88, 0.00) 


5 


(-70.13, 0.00) 


G 


(-72.25, 0.00) 


7 


(-76.44, 0.00) 


8,9 


(-98.22, ±43.87) 


10 


(-106.0, 0.00) 


11 


(-135.4, 0.00) 


12,13 


(-135.9, ±9.318) 



test field method developed by ISchrinner et alj (120051 12007ft 



For a set of suitably chosen test fields, the induction equa- 
tion for the fluctuating field is solved with the actual flow of 
the numerical model, the electromotive forces are computed, 
and then averaged in azimuth and time. Finally the coeffi- 
cients are derived from the projection of the mean electromo- 
ti ve forces on the corresponding test fields. For details we refer 
to Schrinner et al. I d2007l) : ISchrinnerl d201 ll) . For axisymmetric 
mean fields the ct-tensor has 9 non-zero elements while the f3- 
tensorhas 18 non-zero entries. As an example we show in Fig. [3] 
two components for the dynamo model under consideration: 
a w and D w = i] + \fi vr e - \P v e r , i-e. the ^-component 
of the turbulent diffusivity plus the molecular diffusivity 77. For 
more detai ls on the representatio n of the dynamo coefficients 
we refer to lSchrinner et alJ (120071) . in particular Eqs. (10)-(13) 
and (15e). The dynamo coefficients reflect the action of the con- 
vection columns on the mean magnetic field outside the inner 
core tangent cylinder. They are almost zero inside the tangent 
cylinder where no convection takes place. 
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Figure 4: The magnetic field B° (top), the vector potential j4° (middle), and 
the adjoint current J° (bottom) of the fundamental dipole mo de of the numer- 
ical m odel. In the gauge we adopted for the vector potential I Schrinner et al.. 
2010b]), A° is zero due to axisymmetry. 



3.3. Dynamo modes 

We next solve the eigenvalue problem of the dynamo equa- 



tion 



\B = V x DB , (2) 



with 

DB = vxB + aB /3: (VB) - ?yV x B . (3) 

Here v is the mean flow and B the mean field. The eigen- 
functions of the dynamo equation are denoted as B l (r) and 
we refer to them as the 'dynamo modes'. They are constructed 
by expanding them in the magnetic decay modes b (r) of the 
spherical shell V: 



B' 



(4) 



Summation convention over double upper indices is adopted. 
Expansion © is now inserted in Eq. @, after which the eigen- 
values A,; and the matrix elements e li may be determined as 



explained in ISchrinner et al.l (1201 Obi) . The computation of the 
matrix elements requires that the normalisation of B l is fixed, 
and we suppose that that has been taken care of. Details are 



given in Appendix _A Th e advan tage of the new method de- 
scribed by Schrinn er et al. ( 2010bl) is that the velocity and the 



dynamo coefficients do not have to be differentiated, as would 
be necessary in the traditional methods. This is important for 
the present study, as we obtain these coefficients and the veloc- 
ity from the dynamo code in the form of numerical tables. 

Since the dynamo modes B l are in general not orthogonal 
we make use of the adjoint modes B k . Together they form a 
biorthogonal set on V + E: 



V+E 



B k ■ B i d 3 r = I J" -A 1 d 3 r = S k 



(5) 



The second relation follows by integrating by parts, and shows 
that currents and vector potentials also form a biorthogonal set 
on the volume V of the dynamo. The adjoint modes may again 
be found with the help of an expansion in decay modes: 



B = f 



fcs* IS* 



(6) 



and Schrinner et al. ( 20 1 Obh show that f = (e -1 )^ where t in- 
dicates the Hermitean adjoint. For details and the question of 
gauge invariance we refer to this paper. Since the eigenmodes 
and their adjoints are constructed as linear combinations of the 
decay modes they automatically satisfy the boundary condi- 
tions, because the decay modes do so. We refer to Section [XTl 
for a summary of the various magnetic field related concepts 
that we have introduced (model field B, dynamo modes £T, 
decay modes 6\ their adjoints, currents and vector potentials). 

As explained above, we restrict the mode decomposition to 
antisymmetric and axisymmetric dynamo modes. The eigen- 
values of the first of these modes are shown in Table Q] The 
fundamental mode B° is non-periodic and antisymmetric with 
respect to the equator, and may, on that account, be represen- 
tative for the geodynamo. Non-periodic modes B l are them- 
selves real. The structure of the fundamental mode B° together 



4 



with the corresponding vector potential A and adjoint current 
J° is shown in Fig. [4] 

Remarkable is the large drop in growth rate beyond the fun- 
damental mode which underlines the outstanding role of the 
fundamental mode in representing the magnetic field of this dy- 
namo model. For completeness we mention that the spectrum 
of the symmetric modes is similar to that of the antisymmet- 
ric modes, but the fundamental symmetric mode decays more 
strongly, at a rate of —6.298 rj/L 2 . 

3.4. Field decomposition 

The magnetic field B of the numerical dynamo model is rep- 
resented by an expansion in dynamo eigenmodes B l : 

B{r,t)=a\t)B\r) (7) 

So we first construct the dynamo modes from the decay modes, 
and then we represent the field B of the dynamo as a superposi- 
tion of dynamo modes. Because the dynamo modes (eigenfunc- 
tions of the dynamo equation) may be supposed to constitute a 
complete function set they may be used to represent any arbi- 
trary magnetic field. 

The expansion coefficients a 1 may be computed by taking 
the inner product with the adjoint B k and integrating over all 
space: 



B 



Bd 3 i 



V+E 



Ad 3 i 



(8) 



Since B° is real, the coefficient a is real. To evaluate the inte- 
grals in ([8]) we compute first the adjoint currents J k = { ki *j * 
on the numerical grid, cf. Eq. (5). This needs to be done only 
once. Next we compute at each time step the vector potential 
A. The integral f v J k - A d 3 r in ([H) follows by taking the inner 
product and summing over the grid. Because of the infinite vol- 
ume and the power-law scaling of B k ■ B with radius, it is not 
a good idea to evaluate the other expression fy+E B k ■ B d 3 r. 
Fig.|5]shows |a fc (i)| for k — and k — 1 as a function of time, 
together with the averages (\a k \) (horizontal lines). Once the 
expansion coefficients have been determined the computation 
of the cross correlation coefficient (a k a e *) is straightforward. 

Two remarks on the magnitude of the expansion coefficients. 
The a k are of order unity, and that is a result of scaling: in the 
dynamo code B has been scaled to be of order unity. And the 
normalisation ( IA.U makes that the integral in Eq. $Q is also of 
order unity, at least for the modes of lower order. The second 
remark is that by adopting the normalisation ( 1A.1I ) we ensure 
that all dynamo modes, if excited alone and with the same am- 



plitude (i.e. c 
energy, see Appendix A 



1 other a 1 = 0), have the same total magnetic 



4. Theoretical mode excitation levels 



The r.m.s. mode coefficients may be found from the evolu- 
tion e quation for the correlation coefficients derived in Hovng 
d2009h : 
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f ) (a m a n *) . (9) 
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Figure 5: Modulus of the amplitude of the fundamental mode (k = 0, top) 
and of the first overtone (k = 1, broken line), as a function of time in units 
of the magnetic diffusion time L 2 /r]. The horizontal line indicates the mean 
modulus. 
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Figure 6: The cross correlation function {C km (i)C en * (t — r)) as a function 
of the time shift r, for k = I = 1, m = n = 0. The vertical axis is in units of 
WL 2 ) 2 =4(r,/L 2 ) 2 . 



(summation over n, m, not over k, t). The matrix elements 
M kmtn are defined as integrals over velocity cross correlation 
functions: 



M 



km ( i- 



dT{C km (t)C en *{t-T)) . 



(10) 



The functions C km (t) are in turn averages of the turbulent flow 
u over the volume of the dynamo weighted by 'spatial filters' 
defined by the dynamo modes: 



C kt (t) 



u -J k y.B l d 6 r 



(11) 



This expression suggests that C ke may be interpreted as being 
proportional to the work done per unit time by a part of the 
Lorentz force on the flow u. 



4.1. Computation of M kmln 

The turbulent flow u must be measured and projected on 
a set of known vectors J k x B l , and summed over the spa- 
tial grid. This may be done on line after each new time step 
of the dynamo code. The result would be a set of time series 
C ke (t). However, the mean flow v and therefore the turbulent 
flow u = V — v is not known until the end of the simulation, 
so we must adopt a more elaborate procedure. We project the 
full flow field V on J k x B e , and define C krn as above but with 
the full flow field. The required cross correlation may now be 
computed from 

(C km {t)C in *{t-T)) 

= (C km (t)C en *(t - r)> - (C km )(C en *) . 

(12) 

A result is shown in Fig. [6] The computation of the averages 
(C ) in (fT2l is straightforward, and we briefly explain how 
the cross correlation functions are computed, suppressing all 
details on the implementation in the dynamo code. Writing mo- 
mentarily x = C km and y — C n *, we compute, for a given 
time t, an array of values x(t)y(t — i ■ At) for i — 0, ..n. Next 
we select other (later) starting times t, iV av in total, and average 
x(t)y(t — i ■ At) over these starting times t for fixed i. In this 
way we obtain (x(t)y(t — r)) for time shifts r ranging from to 
t cu = nAr. The upper limit r cu to the various correlation times 
t c is not known a priori but has to be determined empirically 
from graphs similar to Fig. [6] The results presented here have 
been derived with r cu = 0.113 i 2 /^ Ar = 7.5 • 10" 4 i 2 / ? 7 
and N av = 6 • 10 4 . Finally we determine M kmin by an inte- 
gration over r, in accordance with Eq. ([Tol l, carried out in the 
post-processing. 

5. Results 

We shall now compare the numerical results with the theo- 
retical predictions and draw some conclusions. 

5.1. Applicability of the statistical theory 

The statistical theory that we used here has been formulated 
by Hovngl ( 2009h . and the conditions of the applicability of the 
theory have been worked out in Sections V and VI of that paper. 
The first requirement is that the correlation time t c is short: 



C U T C « 1 



(13) 



In Hovngl ( 2009 ) C M could be estimated only approximately, 
but here we can do better with the help of the numerical 
model. From Fig. [6] at r = we deduce that C 10 ~ [38-4 
(ry/L 2 ) 2 ] 1 / 2 ~ 12 rj/L 2 . The correlation time may be taken as 
the 1/e time in Fig. [6] r c ~ 0.01 L 2 jr\. In this way we obtain: 



C 



10, 



0.1 



(14) 



C can be shown to increase slowly with mode number k as 
L/Lk, where Lj~ is the characteristic length scale in mode k 

(L k <L = L). 
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Figure 7: Mode excitation levels (\a k \ 2 ) as a function of the mode number 
defined in Table 1: * = (|a fe | 2 ) as measured from the numerical model; = 
( | a k | 2 ) computed from j!6> ; A = ( | a k | 2 ) computed from {5) . The excitation 
levels have been normalised to that of the fundamental mode, whose value is 
measured to be (|a°| 2 ) = 0.226. The levels of mode 1 and 2 are equal since 
they are a complex pair; likewise, the levels of mode 8 and 9 are equal. 



The second requirement is that the unperturbed evolution of 
the system is slow. Since the mean flow v is very small, thi s 
condition reduces to the requirement derived in lHovng (2009), 
relation (21): 



0.01 



L 2 



LkLi 



< 1 



(15) 



From these results we conclude that the first few dynamo modes 
of the numerical model obey the condition of a short correlation 
time and of a slow unperturbed evolution, though not by a wide 
margin. We return to this point in Section fOl 

5.2. Mode excitation levels 

We now come to the r.m.s. mode excitation levels. We con- 
sider first the special case that the fundamental mode is dom- 
inant, that is, (|a°| 2 ) <C (|a 4 | 2 ) for i > 1. In that case we 
may ignore all terms in the summation in the right hand side of 
Eq. ((9]) except m = n = 0. We put d/dt = (the dynamo 
is supposed to be in a quasi-stationary state) and take k = £ to 
obtainQ 



<|aT> 



(16) 



Once the M hoko have been computed as in Section l4"Tl relation 
( TToT l allows us to obtain the overtone excitation levels relative to 
that of the fundamental mode. In much the same way we could 
also derive expressions for the cross correlations (a m a n *), but 
we restrict ourselves here to the autocorrelations. The result is 



'in relation (40) of lHovnd <2009h . the 5R in front of M koko has been erro- 
neously omitted. 
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plotted in Fig. [7] as 0, while the measured overtone excitation 
levels are plotted as *. The agreement between theory and mea- 
surement seems reasonable to good for all modes except mode 
5, 8 and 9. 

As it is not clear to what extent the fundamental mode am- 
plitude is really dominant, we have investigated if the solution 
improves by employing the full equation @. To this end we 
take again k = £ and d/di = in Eq. (0 and evaluate the entire 
right hand side, the M kmln as in Section @j] and the (a m a n *) 
are measured from the numerical data. In the end we divide by 
— 25RAfe. The resulting excitation levels have been plotted as A 
in Fig. [7] Since the A's coincide closely with the 0's, our tenta- 
tive conclusion is that the higher order terms in the summation 
in Eq. (0 may indeed be ignored, and that fundamental mode 
amplitude of the numerical model is in fact dominant. 

The 'flat leve l' in Fig . [7] may be explained with the help of 
estimate (43) of lHovn g (2009): 



(Kf) 

(| a 0|2) 



N~ 



0.1 



(17) 



which agrees rather well with Fig. [7] Here N is the number 
of convective elements, which we estimate to be of order 10 
(Fig. |2). So while there is a good agreement at the qualitative 
level, we have no stringent explanation for the small measured 
excitation levels of modes 5, 8 and 9. However, we recall that 
assumptions (TPfl i and (TT~5T > are only marginally fulfilled, in par- 
ticular for higher overtones. Moreover, we may have lost ac- 
curacy in each of the computational steps visualised in Fig. Q] 
Notably, the parametrisation of the electromotive force by the 
dynamo coefficients determined and thus the dynamo modes 
and their eigenvalues will be to some ex tent inaccurate (see also 
Schrinner et ail 120071: ISchrinnenl201 lh . 

It is tempting to interprete Fig. [7] as a spectrum, but the non- 
orthogonality of the dynamo modes renders such an interpreta- 



tion impossible, see Appendix A for more details 




1.10 
-1.10 



Figure 8: The bottom row shows the r, 8 and tp component of the dynamo field 
B averaged over azimuth and time. The rows above the horizontal line show 
the field expansion S^~q a k B k averaged over the same period and truncated 
at various levels. At the top is p = 0, in the downward direction followed by 
p = 5, p = 7 and p = 13. 



5.3. Representation of the field 

The fundamental dynamo mode largely dominates the mag- 
netic field as illustrated in Fig. [8] It is obvious that in particular 
the r and 9 components of the time and azimuthally averaged 
dynamo field B (bottom row) are fairly well represented by the 
fundamental dynamo mode alone (top row). The ip component 
improves somewhat by including a few overtones, but as we do 
so the r and 9 components deteriorate slightly. This problem 
is fixed again by including a few more modes. By the time we 
have incorporated some 10 modes the agreement is quite good. 

The behaviour of the ip component shows a peculiarity in- 
herent to an expansion in non-orthogonal modes: the approxi- 
mation is not monotonous. Another illustration is the fact that 
mode 6 and 7 have a large amplitude, see Fig. [7] but they almost 
cancel each other in the field expansion ©, and together they 
contribute little to the representation of the field, as can be seen 
in the middle row of Fig. [8] 
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6. Discussion and summary 

Dipolar geodynamo models in the low Rossby number 
regime are governed by one fundamental dynamo mode. On 
average, the fundamental mode has zero growth rate, whereas 
more spatially structured overtones are highly diffusive. Small- 
scale contributions to the magnetic field (i.e. overtones) result 
from the deformation of the fundamental mode by the turbulent 
flow; subsequently, they decay due to Ohmic diffusion. The 
ratio between both processes determines their excitatio n lev- 
els. T his intuitive picture proposed and formalised by Hovngl 
(2009) is confirmed by the satisfactory agreement between the- 
oretically predicted and numerically measured excitation levels 
in Fig. [7] Moreover, from a methodological point of view, the 
consistency of the whole 'train' of computational steps illus- 
trated in Fig.[T] is confirmed. 



As noted by ISchrinner et al.l d2010al) . the dominance of the 
fundamental mode in models of this dynamo regime is closely 
related to their time dependence and to the saturation of the 
magnetic field. The dipole field in these models is stable and 
polarity reversals are inhibited, because the damping of the 
magnetic field results in the preference of always the same 
fundamental mode, which is continuously quenched and re- 
built. As a b y-product, these mode ls may be treated kine- 



matically (see Sch rinner et all |2010a) and the mode analysis 
presented in this study is applicable. However, a generalisa- 
tion of our approach to models beyond the kinamatically sta- 
ble regime would require to place v elocity modes and mag- 



netic field modes on equal footing dCourvoisier et all 120 lot 



Rheinhar dt and Brandenburg , 201C ) and is therefore not easily 
feasible. 

We expand the magnetic field in dynamo modes which are 
eigenfunctions of a non-selfadjoint operator. They form a com- 
plete function set but are in general not orthogonal. This has 
severe consequences for the interpretation of Fig. [5] or Fig. |7J 
they cannot be interpreted as spectra. In general, there is no 
correspondence between the modulus of a mode coefficient and 
the contribution of its associated mode to the dynamo field B. 
Modes with a comparatively large amplitude may in fact con- 
tribute only very little. This has been demonstrated in Fig.[8]for 
modes 6 and 7 and is the reason for the non-monotonous series 
approximation ||7). The ambiguity of mode coefficients is a gen- 
eral problem inherent to most of the attempts to model dynamo 
action by amplitude equations. It could be avoided by expand- 
ing B in a orthogonal function set, e.g. the decay modes, in- 
stead of the dynamo modes. On the other hand, dynamo modes 
lead to a rather accurate representation of B with a minimum 
number of modes taken into account. Furthermore, the com- 
pact formulas obtained for the excitation levels, (0 and ( fT6l l. 
requir e that B i s exp anded in eigenfunctions of the dynamo op- 
erator |Hoyngl [2009). In short, an expansion in dynamo modes 
has been preferred because it leads to the simplest description 
and allows for theoretical predictions. 

In this s t udy, w e have applied the mode analysis developed 
by Hovngl ( 20091) to a geodynamo model in the low Rossby 
number regime. We find satisfactory agreement between the- 
oretical predictions and numerical simulations. The expansion 



of B in dynamo modes may serve as a useful analytical tool 
to interpret numerical dynamo models. Besides direct numer- 
ical simulations, low-order models have been successful in il- 
lustrating the possible interaction of large- scale modes of dif- 
ferent symmetrie s in planetary cores (e.g. lPetrelis et al. , 2009; 



Gissinger, I20101) . The work presented here connects both ap- 



proaches and allows for the possibility to compare predictions 
from low-order models with direct numerical simulations. 
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Appendix A. Mode normalisation and spectral interpreta- 
tion 

First of all, we note that relation (O does not fix the nor- 
malisation of the dynamo modes: given a biorthogonal set 
B k ,B l we may construct others, viz.: B k -> [iB k , B k — > 
(1/ n)B k , where the factor /i may depend on the mode number 
k. These new sets also obey relation Q. We choose to remove 
the remaining ambiguity, up to an overall phase factor exp(i0), 
by imposing that 



B k * ■ B k d 3 r 



V+E 



J k * -A k d 3 r = 1 (A.l) 



(no summation over k). The middle relation follows as usual 
by integrating by parts. An immediate consequence of this nor- 
malisation is that 



(ee t ) fefc = 1 , 



(A.2) 



(no summation over k). The diagonal elements of ee' are all 
unity. To see this insert relation (0]i in eq. ( IA.U and use the 
orthogonality of the decay modes. 

The normalisation ( IA. 1 b makes sense for two reasons. The 
first is that in the special case of self-adjoint modes we may 
choose B k = B k * , and then relation ( IA.U follows automat- 
ically from ©. The second reason is that, loosely speaking, 
dynamo modes with the same amplitude a k but excited alone, 
contain the same total magnetic energy. We shall now make this 
statement precise. 

The energy density of the magnetic field may be written as 
B ■ B = a i a k *B i ■ B k *, according to Eq. ©. Next, we insert 
twice relation ©, integrate over all space, and use the orthogo- 
nality of the decay modes. Finally we apply relation ( IA.2I ). The 
result is: 



f B-Bd 3 r = V 

JV+E u 



,k\2 



£ ( eet )' 



k a i a k * 



(A3) 



This expression shows that if only one mode is excited at unit 
amplitude, a J = 1 and all other a k zero, then the total magnetic 
energy f v+E B ■ £> d 3 r is also unity, independent of the mode 
number j. 
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Now if the dynamo modes were self-adjoint, then Eq. © 
becomes f v+E B k * ■ B i d 3 r = 5 hl . Upon inserting again 
relation we get ee^ = I, i.e. the matrix c is now unitary. In 
that case the second term in (1A.31 > drops out, and we are left with 
fv+E B ■ B d 3 r = J2k \ ak \ 2 - And that leads straight away to a 
spectral interpretation of Fig. [7] because the mode number scale 
on the horizontal axis can be converted into a wave number 
scale. But in general we are stuck with Eq. ( lA.3t as a whole, 
and contributions to the magnetic energy remain stored in the 
nonzero mode overlaps f v+E B k * ■ B l d 3 r. 
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